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Influence of trap location on the efficiency of trapping in dendrimers and 
regular hyperbranched polymers 
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The trapping process in polymer systems constitutes a fundamental mechanism for various other dynamical 
processes taking place in these systems. In this paper, we study the trapping problem in two representative 
polymer networks, Cayley trees and Vicsek fractals, which separately model dendrimers and regular hyper- 
branched polymers. Our goal is to explore the impact of trap location on the efficiency of trapping in these 
two important polymer systems, with the efficiency being measured by the average trapping time (ATT) that 
is the average of source-to-trap mean first-passage time over every staring point in the whole networks. For 
Cayley trees, we derive an exact analytic formula for the ATT to an arbitrary trap node, based on which we 
further obtain the explicit expression of ATT for the case that the trap is uniformly distributed. For Vicsek 
fractals, we provide the closed-form solution for ATT to a peripheral node farthest from the central node, as 
well as the numerical solutions for the case when the trap is placed on other nodes. Moreover, we derive the 
exact formula for the ATT corresponding to the trapping problem when the trap has an uniform distribution 
over all nodes. Our results show that the influence of trap location on the trapping efficiency is completely 
different for the two polymer networks. In Cayley trees, the leading scaling of ATT increases with the shortest 
distance between the trap and the central node, implying that trap's position has an essential impact on the 
trapping efficiency; while in Vicsek fractals, the effect of location of the trap is negligible, since the dominant 
behavior of ATT is identical, respective of the location where the trap is placed. We also present that for all 
cases of trapping problems being studied, the trapping process is more efficient in Cayley trees than in Vicsek 
fractals. We demonstrate that all differences related to trapping in the two polymer systems are rooted in 
their underlying topological structures. 

PACS numbers: 36.20.-r, 05.40.Fb, 05.60.Cd 



I. INTRODUCTION 

As a paradigmatic dynamical process, trapping de- 
scribes and characterizes various other important phys- 
ical processes on complex systems, e.g., page search or 
access in the World Wide Webi"— . The trapping problem 
constitutes an integral primary problem of random walks, 
defined as a kind of isotropic random walks with a perfect 
trap located at a given position, absorbing all particles 
that visit it^. The highly desirable quantity related to the 
trapping problem is the trapping time (TT), also known 
as mean first-passage time (MFPT)^"^'^, which represents 
the expected time for a walker starting off from a source 
node to arrive at the trap for the first time. The average 
of trapping time over all starting nodes is defined as the 
average trapping time (ATT), which is very important 
and useful in related fields, since it is often used as a 
quantitative indicator measuring the efficiency of trap- 
ping process. 

In view of the theoretical and impractical relevance, de- 
termining ATT in diverse systems has received a tremen- 
dous amount of attention within the scientific commu- 
nity. Thus far, trapping problem has been extensively 
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studied for a lot of complex systems, such as regular 
lattices with different dimensions^iiir— , the treelike T— 
fractalsi^"— , the small- world uniform recursive trees^i^^, 
the Sierpinski gasketa^^"— , as well as fractal^^""— or non- 
fractal22r— even modular— scale-free graphs. These 
works have unveiled some nontrivial infiuences of certain 
particular structural properties on the leading behaviors 
of ATT for trapping problem performed in these different 
systems. 

In addition to aforementioned systems, trapping in 
polymer systems'^^ is also a fundamental topic due to 
its wide range of applications, including lighting har- 
vesting in antenna systems^^— , energy or exciton trans- 
port in polymer systems^i^, and so on. Among vari- 
ous polymer systems, Cayley trees^i^ and Vicsek frac- 
tals^"— are two important classes modeling separately 
dendrimers and regular hyperbranched macromolecules, 
both of which have been and continue to be active sub- 
jects of research in numerous fields^^"^®. Recently trap- 
ping problem in Cayley trees^i^i^ and Vicsek fractals^ 
with a deep trap at the central node has been studied, 
uncovering how the underlying structures of the two poly- 
mer networks affect the efficiency of this special trapping 
process on them. However, previous work a^^'^^" — have 
shown that for trapping in a network (e.g., scale- free 
network^SiZS) , the trapping efficiency may rely on the 
position of the trap. Then, two interesting open ques- 
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tions arise naturally for Cayley trees and Vicsek fractals: 
What is the impact of trap's location on the behavior for 
ATT in these two polymer networks? And how the ATT 
scales with the system size when the trap is uniformly 
distributed over the networks? 

In this paper, in order to explore the role of trap's po- 
sition on the trapping efficiency in polymer systems, we 
study the trapping problem on Cayley trees^i^ and Vic- 
sek fractals^"—, both of which grow in an iterative man- 
ner and have a treelike structure that allow for analyti- 
cally treating most cases of trapping processes occurring 
on them, obtaining explicit expressions for the ATT char- 
acterizing the trapping processes. We focus on attacking 
two cases of trapping problems on the two polymer net- 
works. In the first case, the trap is fixed on a particu- 
lar node, while in the other case, the trap is distributed 
uniformly over all nodes in separate networks under con- 
sideration. We obtain analytical closed-form solutions or 
numerical solutions to the ATT for both cases. 

Concretely, for Cayley trees, we derive an exact for- 
mula for ATT when the trap is located at an arbitrary 
node, the leading behavior of which depends on the short- 
est distance between the trap and the central node: the 
longer the distance, the higher the ATT, indicating that 
the position of trap has a substantial effect on the trap- 
ping efficiency. While for Vicsek fractals, we obtain an 
exact analytical solution to ATT when the trap is a pe- 
ripheral node, as well the numerical solutions to ATT 
when the trap is placed on another node, with both so- 
lutions having the same dominant scaling, implying that 
the influence of trap's location on the ATT is negligible. 
For the second trapping problem, by using different ap- 
proaches, we deduce accurate expressions for the ATT 
for both networks, which show the trapping efficiency of 
Cayley trees is much higher than that of Vicsek fractals. 
We present that their particular structures are responsi- 
ble for the difference of results obtained for the two net- 
works. This work deepens the understanding of trapping 
problems in the two representative polymer networks. 



II. NETWORK CONSTRUCTIONS AND PROPERTIES 

In this section, we introduce the two representative 
polymer networks — Cayley trees and Vicsek fractals — 
and their structural properties. Both networks are con- 
structed in a deterministically iterative manner. 




FIG. 1. (Color online) Structure of the a particular Cayley 
tree 6*3,5 • 



we add to — 1 new nodes and link them to the peripheral 
node. All the new introduced nodes at this stage become 
the peripheral nodes of Cm,g- Figure [1] illustrates the 
construction process for a specific Cayley tree C3.5. 

By construction, one can easily check that at each it- 
eration i {i > 1) the number of newly created nodes is 
Ni{g) — m{m — 1)'^^. Thus, for Cm,g, the number of 
peripheral nodes and the total number of nodes are 



(1) 



and 



iVg = l + £iV,(g) = 



to(to — ly 



(2) 



respectively. Then, the total number of edges in C'm,g 

m[m — 1)3 — TO 



IS 



N„ 



1 



(3) 



It should mentioned that although Cayley trees ex- 
hibit an obvious self-similar structure, they are nonfrac- 
tal since they have an infinite fractal dimension. 



B. Vicsek fractals 



A. Cayley trees 

Cayley trees^^^ after g iterations, denoted by C,n,g 
(to > 3,3 > 0), are built in the following iterative way. 
For 5 = 0, Cm,o comprises only a central node. For g = 1, 
TO nodes are generated connecting the central node to 
form Cm,i, with the to single-degree nodes constituting 
the peripheral nodes of Cm,i- For any g > 1, Cm,g is ob- 
tained from Crn.g-i- for each peripheral node of Cm.g-i, 



Vicsek fractals^"— after g iterations, denoted by Vf^g 
if > 3, g > 0), are constructed in a different iterative way 
from that of Cayley trees. For g = 0, V/,o consists of an 
isolated node without any edge. For g = 1, f new nodes 
are generated with each being connected to the node in 
Vffl to form Vf^i, which is exactly a star. For g >2, Vf^g 
is obtained from Vjg_i. To obtain Vf^g, we introduce / 
new identical copies of Vf^g-i and arrange them around 
the periphery of the original Vf^g-i. Then we add / 
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FIG. 2. Illustration of the first several iterative processes of 
a particular Vicsek fractal, V4,3. 




FIG. 3. (Color online) Another construction method of Vicske 
fractals, u and v are two adjacent nodes in Vf^g-i. At iter- 
ation g, each of them gives rise to / new neighbors, denoted 
by u\, U2,. ■ Uf, and «i, V2,. ■ Vf, respectively. Since u 
and V are directly connected by an edge in V/.g-i, two of 
their new neighboring nodes (e.g., ui and Vi) are connected 
to each other by a new edge. 



new edges, each of them connecting a peripheral node 
in one of the / corner copy structures and a peripheral 
node of the original central structure, where a peripheral 
node is a node farthest from the central node. Figure [2] 
shows the structure of a special Vicsek fractal 14^3. By 
construction, at each generation the number of the nodes 
increases by a factor of / + 1; therefore, the total number 
of nodes of Vf^g is Ng = {f + 1)^, and the total number of 
edges in V/,g is £'g = A'^g — 1 = (/ + 1)^ — 1. In contrast to 
the Cayley trees, Vicsek fractals are fractal objects with 
the fractal dimension being equal to log3(/ + 1). 

We continue to study some properties of the Vicsek 
fractals, which are useful for the following text. Let Lg 
denote the diameter of Vf,g . It is easy to see that Lg sat- 
isfies recursive relation Lg = 3Lg-i + 2, which together 
with the initial condition Li ~ 2 leads to Lg — 3^ — 1. Ac- 
tually, for two nodes in Vf^g, if their shortest distance is 
equal to the diameter Lg , then the shortest path connect- 
ing the two nodes must include the central node of Vf^g , 
and these two nodes are two peripheral nodes of Vf^g . Let 
Ng denote the number of peripheral nodes of V/,g. Evi- 
dently, Ng satisfies recursive relation Ng = (/ — l)iVg_i. 
Considering Ni = f, we have 

Ng = fif - iy~' . (4) 

In addition to the above replication and connection 
operations, Vicsek fractals can also be alternatively con- 
structed using another method^^ as illustrated in Fig. [31 
Suppose one has Vf^g-i. Then, Vf^g can be obtained from 
Vf^g^i by performing the following steps. First, for each 
node in V/_g_i, / new nodes are generated and linked to 
the old node. Then, for each pair of adjacent nodes, u 
and V in Vf^g-i^ a new edge is added between two of their 
new neighboring nodes. Moreover, each new node gener- 



ated at generation g can have at most one new neighbor 
born at the same generation. 

According to the second construction way, we can cate- 
gorize the nodes in Vf^g. Let Ag represent the set of nodes 
of Vf^g, and Ag the set of those nodes of Vf^g, which are 
created at iteration 5. Obviously, Ag = Ag_i+Ag. More- 
over, we can classify the set Ag into two subsets Ag^'' and 
A^^\ such that Ag = A^^^ U A^^\ where A^^^ is the set of 
nodes with degree 1 and Ag is that of nodes with degree 
2. It is easy to verify that the numbers of nodes in Ag^"* 
and Ag are 

|A«l = (/-2)(/ + ir^+2 (5) 

and 

|Af)H2(/ + l)s-i-2, (6) 

respectively. 

We now study a new quantity Dg, which represents the 
distance to a given peripheral node in Vf g and is defined 
by 

Dg= ^^(5), (7) 

where di {g) is the length of the unique shortest path from 
node i to the peripheral node in V/_g. According to the 
self-similar structure of V/_g, we have the following recur- 
sive relation: 

Dg = Dg^l + {Dg^l + Ng^l + iVg_lLg_l) 

+ {f - l){Dg-i + 2Ng^i + 2iVg_lLg_l). (8) 
Considering Z3i = 2/ — 1 , Eq. ([5]) is solved to yield 

Dg^\{i<^-mf-l){f+iy-\ (9) 
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After introducing the two polymer networks, in what 
follows we will study the trapping problem defined on 
them, with an aim to gain a comprehensive understand- 
ing of the effect of trap's location on the absorbing effi- 
ciency on Cm,g and Vf^g. For this purpose, we will first 
investigate random walks with a single immobile trap 
fixed on a certain node; then we will continue to address 
random walks with the trap distributed uniformly over 
all nodes. 



III. TRAPPING WITH AN IMMOBILE TRAP 

In this section, we study isotropic random walks with 
a single trap defined on Cm.g and V/_g, respectively. The 
random-walk model considering here is a simple one. At 
each discrete time step, the walker moves from its current 
location to an arbitrary nearest neighbors with the same 
probability. Let Tij{g) denote the MFPT from node i to 
J, which is the expected time taken by a walker starting 
from i to first reach j. Let Tj{g) denote the ATT to 
trap node j. Then, the interesting quantity related to 
trapping problem is given by 



In the sequel, we will study analytically Tj{g) for both 
Cm.o and Vf,g, and show how Tj[g) varies with network 



size. 



A. Trapping In Cayley trees with a trap at an arbitrary 
node 



In a recent paper—, we have studied the trapping pro- 
cess on Cayley trees with an immobile trap located on 
the central node. Here, we attack a general trapping 
problem on Cayley trees with the trap fixed on an ar- 
bitrary node. To attain this goal, we first classify the 
nodes in Cm.g into g + \ levels according to their shortest 
distance to the central node: The central node is at level 
0, the nodes created at generation 1 are at level 1, and 
so on. For convenience of the following description, let 
ri represent a node at level i, where the subscript i de- 
notes simultaneously the generation it is created at, and 
its shortest distance to the central node. Note that all 
nodes at the same level are equivalent to each other, in 
the sense that the ATT is the same, if any of them is 
considered as a trap. 

Let T^""™(ff) be the sum of the MFPT from a starting 
point to a target node at level j (0 < j < g) in Cm.gi 
where the sum is taken over all starting nodes in Cm.g- 
That is, 



(11) 



Then, the ATT to an arbitrary node at level j in Cm,g is 

1 



TrAg) 



N, 



.T--(g). 



(12) 



Thus, to find Tj. {g), we will alternatively evaluate the 
quantity T^^^'ig). 

In order to obtain T™"^((7), we regard Cm,g as a rooted 
tree with the central node being the root. Then, the 
following relation holds: 



Tc"'i9) + - — Tr,,r, (g) 



m 



771 — 2 



m - 2 



Tr\g-^) + im-ir-^Tr^rAg) 



777 — 1 

H Tc [g-]), 



(13) 



where Tc""'(5) the sum of MFPTs for all nodes to the 
central node on Cm,g, and Tr-r^ is the MFPT from a node 
at level i to one of its offspring node at level j. 

Equation (fT3|) can be accounted for as follows. The 
first term is based on the fact that a walker starting from 
a node, which and the trap node have the only one low- 
est common ancestor (i.e., the central node), should first 
(10) visit to the central node, and then takes Trorjig) more 
time steps to reach the trap. Here, the lowest common 
ancestor for two nodes i and j is defined as the node 
with the possible biggest level value in the rooted tree 
but having both i and j as its descendants. The second 
term describes the case that a particle starting off from a 
node, which and the trap have the lowest common ances- 
tor at level i, first jumps to the lowest common ancestor 
and then takes Tr-r- (<?) more steps to arrive at the target 
for the first time. The last term accounts for the sum 
of MFPTs from all descendants of the trap to the trap 
itself. 

Next we derive the two quantities T™™((7) and Tr-r, (g) 
with i < j. In a previous work^, we have derived that 
the MFPT, Tr^cig), from a node at level i to the central 
node is 



Tr,c(5) 



(m - 2y 



g-t+ll 



(14) 



utilizing which we can easily determine T™™^{g) given by 



r6"'"(5) = £^^(3)r.,c(.9) 



(to - 2)3 



[2(to- 1)29+1 



m{rn + 2)g + m 

(to - 2)2 ^ ^ 



(15) 



We proceed to evaluate Tr.r-{g) {i < j), which can be 
expressed as 



3-1 

Tr,r,{g) =^Tr,r^^,{g). 

k—i 



(16) 
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Before determining Tr^r, (g), we first calculate the MFPT, 
T^-n^iig), from a node at level i (0 < i < g) to its 
neighboring node at level i + 1. For the case of i = 0, it 
is easy to have 



is reduced to the previous result^. For another limiting 
case {j = g) that a specific peripheral node is looked 
upon the trap, by substituting j = g into Eq. (f23| . we 
obtain the ATT as 



1 

m 



2m^{m - 2)g - 2m(m — 1) 



Substituting Eq. (HH) into Eq. dTT]), we can solve Eq. 
to yield 



Trori (g) = 



1 



[2(m- 



TO - 2 

For i > 1, the following relation holds 
1 1 



(18) 



(rn — 2)2[to(to — 


1)9- 


-2] 


2m{m? — 4)g — 


Sto^ 


+ 4 


(to — 2)2[to(to - 


-ly 


-2] 


3m? + 2m 


-4 




(to — 2)2[m(TO - 


-1)9 


-2] 



(to - 1) 



2s 



1)^ 



nn+i(5) = ^ 1 

mm 

m - 2 



(3) + 7;...+,(g)] 



TO 



l + T;ic(5-j)+T;,r,+i(g)] 



(24) 



We continue to express Tr {g) in terms of the network 
size Ng, in order to obtain the relation governing the two 
quantities. Recalling Eq. ([2]), we have (to — 1)^ = [(to — 
2)Ng + 2\/m and g = [ln((TO-2)7Vg-h2)-lnTO]/ In(m-l). 
These relations allow to recast T^- (5) as a function of Ng-. 



(19) 



Considering the initial condition given by Eq. ([T5|. 
Eq. p9)) can be solved inductively: 



Tr,{g) 



2[jm{m — 2) — TO -|- 1] 

to(to — 2)^Ng 

4 



^r.r.+i ig) 



1 



[2m{m - ly 



2(to- 1)9-* -m] . 

(20) 



to(to - 1)J~i(to - 2)^Ng 
(to^ - 4) [1 - (5 + j)m] + Am 



[{m-2)Ng + 2f 

[{m~2)Ng + 2f 



Inserting Eq. (l20l) into Eq. p6)) yields 

m{j - i) 



2{m 



m{m - 2yNg 
2)[(TO-2)7Vg + 2] 



[(to - 2)Ng + 2] 
,(25) 



T^rAg) 



[2(r 



to(to - 1)^-1 (to - 2)3iVg (m-2)3Arg 

Thus, for a very large system, i.e., Ng 00, we have the 
following expression for the leading term of Tr- (<?): 



(to - 2)2 



[(m - 1)^' 



Plugging Eq. (HSj) and Eq. (12Tj) into Eq 

2TO(m - 1)29 



(21) 
we have 



Tr^ig) + I) Ng. 



(26) 



(to - 2)3 
(m - 1)9 



[jto(to — 2) — to -|- 1 + 2(to — 1)""^ 



(to -2)3 
2(to, + 2)(to 



[{g + j){4m - to3) +m^ + Am- a] 



(m - 2)3 



(m-2)3' 



(22) 



Substituting the expression of Eq. ([22)) into Eq. (IT2|) . we 
obtain the rigorous expression for the ATT to the trap 
node at level j on the gih generation of the Cayley trees: 



2m{m - l)29[jm(TO - 2) - m + I + 2(m - 1)^^^] 
(to - 2)2[to(to, - 2)9 - 2] 
(m - lY[{g + i){Am — m^) -I- to2 -|- 4to - 4] 



(to - 2)2[to(to - 2)9 - 2] 
2(to, + 2) (to, - l)9-i+i + rn2 



(m-2)2[TO(TO-2)9-2] 



(23) 



Equation ([25)) provides an explicit formula for ATT to 
an arbitrary node of Cm,g- For the particular case of 
J = that the trap is fixed on the central node, Eq. ([25)) 



Equation (12^ shows that leading behavior of T,. (17) 
is dependent on the level of trap position, i.e., the dis- 
tance from the central node to the trap. Particularly, 
when j — 0, namely the central node is trap, we have 
Tr-{g) iV^; for j — g that trap is located at a periph- 
eral node, Eq. ([M]) implies Tr^{g) ^ Ng In Ng. Thus, in 
the limit of the large network size Ng, for the trap lo- 
cated on the central node, the ATT grows linearly with 
the increasing network size. However, for the trap fixed 
on a particular peripheral node, the leading asymptotic 
NglnNg dependence of Tr^i^g) with the network size is 
in strong contrast with the linear scaling of Tro (5) with 
N„. 



B. Trapping in Vicsek fractals with the trap at a 
peripheral node 

Different from the case of Cayley trees, for trapping 
problem in Vicsek fractals V/.g, it is very difficult and 
even impossible to determine the exact expression for 
ATT when the trap is located at an arbitrary node. But 
for some cases that the trap is fixed at a particular node, 
the problem can be solved analytically. In a previous 
work^, we have obtained the ATT to the central node 
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on Vf,g. Below we will utilize a similar but a little dif- 
ferent technique to address the trapping process on Vf^g 
with an immobile trap positioned on a specific periph- 
eral node. We will show that, for both cases, the leading 
behavior for ATT is identical. 

Let T,p{g) denote the TT for node i, i.e., the MFPT 
from node i to the trap node. Then, for this case the 
ATT Tp{g) is given by 



Tpig) 



(27) 



In order to determine Tip{g), we introduce two new quan- 
tities for n < g: 



and 



Then, we have 



n K9) 



ieA„ 



TtY>{g) 



TiY>{g)- 



sum 

a 



(g) 



and 



Tisum 
9 



(g) 



rpsum 
9 



(5). 



(28) 



(29) 



(30) 



(31) 



In this way, the problem of determining Tp{g) is reduced 
to finding T™f (5) and T^'"^{g). 

We first deduce the reclusive relation for Tg^f{g). Us- 
ing a similar process as the case that the trap is fixed 
on the central node^, we can derive the following law 
governing the evolution for TT of node i: 

T^p{g) = Hf + l)r,p(.9 - 1) + 3M(g - 1) . (32) 

Note that according to the second construction approach, 
for a particle performing random walks in Vf^g, before 
arriving at the peripheral node as the trap, it must first 
visit its unique neighbor that is in fact a peripheral node 
on V/.g-i. Then, we have 



E 
E 

ieA„_ 



Ttpig) 



[3(/ + l)T,p(.g-l) + 3M(g-l)] 



+ Ng^l[2{Ng 
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1]- 



(33) 



The sum term on the right-hand side (rhs) of Eq. ([55)1 
stands for the time spent by a walker to reach the neigh- 
bor of the trap, while the second term accounts for the 
time steps from the trap's neighbor to the trap. Equa- 
tion ([55)) can be readily simplified to 



For T^"™((7), it can be evaluated as follows. By defini- 
tion. 



Trig) - E T^F{g)+ E T^pi9) 

isAi,' 



(35) 



Applying the approach ii*^, we can evaluate the two 
summation terms on the rhs of Eq. psp . and further 
obtain the following recursive relation for T'g"™(g): 

fr\9) - fT^ig) + 4i?s-i + |AW I ~2Ng + 2 (36) 
Plugging Eqs. ([Ml) and Eq dM]) into Eq. ^ gives 

Tr\9) - 3(/ + l)^r™T(g - 1) + 3/(/ + l)Dg,i 
+A^3_i(/ + 1)(27V,-3) + |AW| 



9-1 



2N„ 



(37) 

Considering T™™(1) = 2/^ — 1 and combining the above- 
obtained results, Eq. ([57| can be solved to obtain 



rpsum 
9 



(.9) 



3/2 + 3/ - 2 



3/24 
2/ 



5/ 

- 1 



2/ 



2 

3^/ 



3S(/ + 1)29 _ (J + 1)29 

6/2 + 5/ + 6 



1)^ 



6/2 + 10/ + 4 



(/ 



(38) 



Inserting Eq. (|38|) into Eq. pop , we arrive at the exact 
formula for the ATT with the trap located at a specified 
peripheral node on V/,g: 



Tp{g) 



3/2 3/ - 2 

3/2 5/ + 2 

2f - 1 
J .39 • 



i'if + iy-if + iy 

6/^ + 5/-f 6 



2f + 2 6/2 10/ 4 



(39) 



which can be further represented as a function of network 



size Ng as 



Tp{g) 



3/2 + 5/ + 2^ ' 



2/-1 
'2/ + 2 



(jVg)l°g3(/+l) 



6/2 + 5/ + 6 
6/2 + 10/ + 4 



.(40) 



From this succinct dependence relation of Tp[g) on net- 
work size Ng, we can find that for large networks, i.e., 
Ng ^ CO, the leading term is 



Tpig) 



(7V^)l + log3(/+l) 



(41) 



which is identical to the behavior of ATT Tc{g) when 
the trap is located at the central node^. 



IV. TRAPPING WITH THE TRAP UNIFORMLY 
DISTRIBUTED 

In Sec. mil we have studied the trapping problem on 



Tg-"iig) - Hf + l)r™'°(.g - 1) + ifDg^i + Ng^i{2Ng - 3) .(j^^ an(i y^^^^ an immobile trap located at a given 

(34)node. In this section, we will study the trapping issue 
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on the two networks with the trap uniformly distributed 
throughout all nodes. 

In this case, what we are concerned with is the quantity 
Tg defined as the average of MFPTs over all pairs of nodes 
in the networks: 



N„ N„ 



(42) 



For convenience, we use Ttotig) to denote the summation 
term on the rhs of Eq. (|^^ : 



Then, 



Ttotig) ^^^T,j{g) . 
i=i j=i 



(43) 



(44) 



which is called global average tapping time (GATT). 

By definition, the quantity GATT involves a double 
average: The first one is over all the starting nodes to a 
given trap node, the second one is the average of the first 
one with the trap having a uniform distribution among 
all nodes. In the sequel, we will analytically study Tg for 
Cayley trees Cm,g and Vicsek fractals V/^g, respectively. 

I 



For convenience, hereafter we also call Tg the ATT in the 
case without confusion. 



A. Cayley trees 

For Cayley trees, we can easily determine Ttotig) and 
Tg by using the intermediary results obtained in Sec- 
tion [III^ For rtot(<7), it obeys the following relation: 



Ttotig) = + J2N,{g)T^^"\g) . 

Substituting Eq. (g^]) into Eq. (gS]), we obtain 
2m2 



(45) 



Ttotig) 



■[(m- l)f - l][2m(m- 1)9 - 1 



(to -2)4' 
+ ((to^ - 2m).g - 2to + l)(m - 1)^^'] 



(46) 



Then, the explicit expression for Tg in Cm.g is 
2to2 



T 



[(m- l)f - 1][2to(to- 1)9 



(TO-2)2[TO(m-l)9-2]2 

+((m^ - 2m)g - 2m + l)(m - l)^^ - 1] , 



(47) 



which can be rewritten in terms of the network size iV^ 
as 



^ 2(A^,-l)[(TO-2)jV, + 2]- 

" TO(TO-2)3(iVg)2 



ln(miVg - 2Ng + 2) - In to 



ln(m - 1) 



(m^ - 2to) - 2to + 1 



4m(iV<, - 1) GmiNg - 1) 



{m-2)^Ng (TO-2)3(iVg)2' 

(48) 



In the limit of infinite network size (i.e., Ng — )■ oo), we 
have the dominating term 

Tg^NglnNg. (49) 

This leading asymptotic Ng In Ng dependence of Tg on 
the network size is equivalent to that of T^^ (5) for trap- 
ping in Cm,g with a peripheral node being a trap, but is 
in marked contrast with the linear scaling of T^q (g) for 
the trapping problem when the central node is the trap. 



B. Vicsek fractals 

For Vicsek fractals, the above method for computing 
Tg in Cm,g is not applicable to that in Vf^g. We next re- 
sort to another method for determining Tg in Vf^g, by us- 
ing the connection^i^^ between resistance distance, also 
refereed to as effective resistance, and MFPTs for ran- 
dom walks on a connected graph. To this end, we view 
Vf,g as an electrical network^^ by considering each edge 



in Vf^g to be a unit resistor—. Let Rijig) be the effec- 
tive resistance between two nodes i and j in the electrical 
network corresponding to V/,g. Then, we have following 
exact relatioii^iiZ^ 

T,,[g)+T,,[g)^2EgR,,ig). (50) 

Using this obtained relation governing MFPTs and effec- 
tive resistance, Eq. p3)) can be recast as 

Ttotig) ^EgY,Y.^-J^9)- (51) 

i=i 3=1 

Equation ([5T|) tells us that if we have a method to 
determine the effective resistance, then we can find the 
quantity Tg. Since Vicsek fractals have a treelike struc- 
ture, the effective resistance Rijig) is exactly the usual 
shortest-path distance between node i and j in V/,g, 
which we denote as Sijig). Then, 

Ng Ng 

Ftot (.g) = E E (ff ) = ^9 Stot ig), (52) 
i=i 3=1 
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where 



Stotig) 



(53) 



i=i j=i 



is actually the Winner index— of Vf^g. 

We continue by showing the procedure of determining 
the total shortest-path distance Stotig), which just equals 
the number of edges in the shortest paths between all 
pairs of nodes in Vf^g. Instead of counting the edges in 
the paths, here we count the paths passing through a 
given edge, and then sum the results of all edges in Vf^g . 
Let be an edge in Vf^g connecting nodes i and j, 

and Eij (g) the number of the shortest paths of different 
node pairs, which pass through (i, j). Let Ni^j(g) be the 
number of nodes in V/,g lying closer to node i than to 
node j, including i itself. Then, 

Stotig) = E E 2N^<Jig)NJ<^ig) 

= 2 E N,^,ig)[Ng-N,^,ig)], (54) 

where Nj^iig) ~ Ng ~ Ni^^jig) was made use of. 

We now apply the relation in Eq. (|54p to deduce 
Stotig)- -f'or this purpose, we classify the edges in Vf,g 
into two sets in the following way. Recalling the second 
construction of Vf,g (see Fig.[3l), an arbitrary edge (m, v) 
in Vf^g-i is replaced by three new edges (u,ui), (uijWi) 

and ivi,v). Let e^^^ denote the set of those edges in V/_g, 
with both endpoints of each edge connecting two nodes 

(2) 

having a degree more than one. And let Cg denote set 
of the remaining edges, all of which have exactly an end- 
point with a single degree. It is easy to derive that the 
numbers of edges in these two sets are 

|eW|=3i?,_i=3(/ + l)^-i-3 (55) 



and 



Making use of Eqs. ((58)) and ([59)) . the quantity Stotig) 
can be represented recursively as 

Stot ig)= E E ^■'^■(5) 

= 3(/ + l)'Stotig - 1) - 4i?s-i + 2|ef )|(7Vg - 1). 

(60) 

Considering S'tot(l) = 2p, Eq. (pD|) can be solved induc- 
tively to obtain 

^tot(5) = ^^T^^K/ + 1)^(3/' - 2/)3^' - 3f + 4/ + 4) 



3/ + 2 
-2(/ + 2)]. 

Then, the analytical expression for Tg in Vf^g is 



(61) 



T, 



EgStotig) _ [(/ + 1)^-1] 



iNg)^ (3/-H2)(/ + l) 



I)tTT[-2(/ + 2) 



+ (./ + 1)^(3/' - 2/)3f - 3/^ + 4/ + 4)] , (62) 
which can be rewritten explicitly in terms of network size 



Ng as 



T, = J^p;^m^ 2/)(A^,)i+'°S3(/+i) 
^7V3(3f + 4/ + 4)-2(/ + 2)]. 



(63) 



When Ng — )■ oo, we have the dominating term for Tg in 



Tg - (Ar^)i+ios3(/+i) ^ 



(64) 



growing as a power-law function of network size Ng, a 
behavior similar to that of Tp(g). 

Note that the above method and process for computing 
Winner index and ATT is general for all trees. We have 
used this approach to calculate Tg for Cayley trees and 
recovered the result in Eq. (|47|) . 



\<i^Q^\ = Eg-\e^P\ = iI -2)if + l)<^-'^ + 2, (56) V. RESULT COMPARISON AND ANALYSIS 



respectively. 

Figure [3] implies that for an edge iu,v) in Vf^g-i, we 
have the following relation for the three edges (m,ui), 
(mi,wi), and (wi,w), in Vf^g-. 

EuuAg) + Eu.vAg) + E-oAg) = 3(/ -i- i)'^E^,ig - 1) - 4. 

(57) 

Summarizing the terms on the rhs of Eq. (j57p over all 
the Eg-i edges in Vf^g-i yields 

E E.jig) = 3(/ + l)'5tot(5 - 1) - ^Eg^i- (58) 



ii,])ee 



(2) 

On the other hand, for edges in , we have 



E i?.,(<?) = 2|ef |(iV,-l). 



(59) 



(ij)e£ 



From the result provided by Eq. (|26p . we can easily 
see that the leading behaviors of Tr^ ig) for different j are 
evidently different. This distinction shows that for trap- 
ping in Cayley trees Cm,g with a deep trap, the trap's 
location has a significant effect on the trapping efhciency 
measured by ATT. Specifically, the dominating scaling 
of ATT grows with the distance from trap to the central 
node: the smaller the distance, the more efficient the 
trapping process. For the case that the trap is the cen- 
tral node, the trapping process is the most efhcient, with 
the scaling of ATT growing linearly with system size Ng ; 
while for the case when the trap is placed at a peripheral 
node, the trapping process is the least efficient, with the 
ATT changing with network size Ng as Ng In Ng . 

In addition, Eq. shows that when the trap is uni- 
formly distributed over the whole Cayley trees, the dom- 
inating scaling of Tg behaves with the network size Ng 
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as N g In Ng. This indicates that the linear scahng of the 
ATT Trg{g) to the central node is not representative of 
the Cayley trees, in the sense that Tg scales larger than 
linearly. In contrast, the scaling of ATT to a peripheral 
node is a representative property for trapping process 
taking place in Cayley trees. 

The different scalings for Tr- (g) lie in the peculiar 
structure of Cayley trees. Figure [1] shows that a Cay- 
ley tree actually consists of m branches (regions) , each of 
which is a subtree with a node at level 1 being its root. 
For j = corresponding to the case when the trap node 
is the central node, the walker, irrespective of its starting 
point, will visit at most one branch before being trapped. 
On the contrary, for j > corresponding the case that 
the trap is located on a node with distance j to the central 
node, the particle, starting from a large group of nodes, 
must first visit the central node and then proceeds from 
the central node along the path ro — ri — ■ ■ ■ — rj-i — rj 
until it is absorbed by the trap. According to previous 
result, the MFPT from a node at level j — 1 to its 
direct neighbor at level j is Tr-_^rj = 2AVj_i<rj — 1, 
which increases with j. Thus, the scaling of Tr (g) grows 
with j, as shown in Eq. (|26p. For instance, for the two 
limiting cases of j = and j = g, Tr„{g) ~ Ng and 
T^r {g) ~ -^glu^g since g ^ InA^g. Finally, for the case 
that the trap is uniformly distributed, the phenomenon 
that the leading term of Tg exhibits the same scaling as 
that of Tr^ {g) can be heuristically understood as follows. 
From Eqs. ^ and ([2]), it is easy to see that in Cm,g the 
fraction of peripheral nodes is about ^Eji thus Tr^lg) 
alone can determine the behavior of Tg . 

Quite different from those in Cayley trees, the domi- 
nating scaling of Tc{g), Tp{g), and Tg in Vicsek fractals 
are identical, all of which scale with Ng as (iVg)^+'°S3(/+i) 
as shown in Eqs. (|4T|) and (|64l) . Moreover, extensive nu- 
merical results also show that for case that the trap is 
placed at other node, the leading behavior for ATT is 
also Ng as {Ngy+^°s^'^f+^K In Fig. HI we report the nu- 
merical results for trapping in Va.g (1 < <? < 7) with the 
trap at different nodes, which shows that the ATTs dis- 
play the same scaling, independent of the trap's location. 
Therefore, the location of the trap has no qualitative ef- 
fect on the scaling of ATT for trapping in Vicsek fractals, 
which is in marked contrast to the trapping process oc- 
curring on Cayley trees. 

The root of identical scaling for ATTs corresponding to 
different positions of the trap is attributed to the struc- 
ture of Vicsek fractals. For the case that the trap is at the 
central node^S, the ATT is the lowest, the leading scaling 
of which varies with network size Ng as ^j^"'"'°S3(/+i)^ Pq^. 
the case that the trap is at another node different from 
the central node, the analysis is similar to that of Cayley 
trees. In this case, to find the target in Vf^g, the walker 
must first visit the central node, taking (iVg)^+'°S3(/+i) 
time steps, then continues to jump towards the trap 
along the unique shortest path from the central node 
to the trap. During the process of jumping along the 
path, the MFPT from a node x to its direct neighbor 
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FIG. 4. (Color online) Average trapping time versus the dis- 
tance from the trap to central node in Vicsek fractals corre- 
sponding to / = 3. 



y is 2Nx<y — 1, which is proportional to Ng. Partic- 
ularly, for the special case when the trap is positioned 
at a peripheral node, the ATT is the highest, which is 
approximately equal to Ng x Lg/2 ~ (Arg)i+'°S3(/+i)^ 
where Lg = 3^? - 1 - (Vg)^°S3(/+i) is the diameter of 
Vf,g, as defined above. Since when the trap is fixed at 
another node other than the central node and periph- 
eral nodes, the ATT is between Tc{g) and Tp{g), and 
thus display the same scaling as Tc{g) and Tp{g). Anal- 
ogously, we can account for the leading asymptotic de- 
pendence (7Vg)^"*''°S3(/+i) of Tg corresponding to the case 
when the trap is uniformly distributed over all nodes. 

From the aforementioned results, it can be seen that 
the trapping efficiency in Cayley trees exhibits rich scal- 
ings in the context of the ATT when the trap is placed 
at different positions. However, the trapping efficiency 
is identical for Vicsek fractals, despite of the trap's loca- 
tion. In addition to this distinction, there are some other 
differences between the ATT for these two networks. For 
example, when the trap is the central node, the trap- 
ping efficiency for Cayley trees is higher than that of 
Vicsek fractals''''. Moreover, by comparing Eq. p6| and 
Eq. (HH), we know that for the case of trap being lo- 
cated at a peripheral node, the ATT of Cayley trees is 
also much smaller than that of Vicsek fractals. Finally, 
for the case that the trap is uniformly distributed, the 
trapping process in Cayley trees is also more efficient in 
Vicsek fractals, which can be seen by comparing Eq. 
with Eq. (|64p . Thus, in comparison with Vicsek fractals, 
Cayley trees have a desirable structure favorable to the 
trapping process. 

It has been proven2S that for trapping problem in a 
general connected graph with an absorbing node, the 
possible minimal scaling for ATT is proportional to the 
graph size and the inverse degree of absorbing node, 
which provides a maximal scaling for the lower bound 
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of ATT for trapping in an arbitrary network with a per- 
fect trap. In this sense, for trapping in Cayley trees, the 
possible minimal scaling of the ATT can only be reached 
when the trap is located at the central node, since in this 
case the leading scaling of the ATT grows proportion- 
ally to the system size and the reciprocal of the degree 
of the central node; while for the case when any other 
node is the trap, the magnitude of the ATT is greater 
than the corresponding possible minimal scaling. This 
phenomenon is in sharp contrast to that observed for 
the hierarchical scale-free grapb^"— , where the possible 
minimal scaling of ATT can be achieved for any node. 
In the context of Vicsek fractals, the possible minimal 
scaling cannot be reached for any trapping node, which 
is also different from that for Cayley trees. 

The above result implies that for trapping in a graph, 
although the possible minimal scaling can be reached for 
a given trap, it may not be achieved when the trap is 
placed on another node. In fact, the possible minimal 
scaling for ATT to a trap depends not only on the posi- 
tion of the trap, but also on the whole topological struc- 
ture of the graphic. In future study, it is interesting to 
explore the problem of designing networks with desirable 
architecture, for trapping process taking place on which 
the possible minimal scaling can be achieved for any node 
as a trap. 



VI. CONCLUSIONS 

We have performed an in-depth study of the trapping 
problem on two classical polymer networks — Cayley trees 
and Vicsek fractals — with a goal to reveal the influence 
of trap's location on the trapping efficiency. To this aim, 
for both networks, we first studied the trapping prob- 
lem with the trap located at a given node; then we ad- 
dressed the case with the trap uniformly distributed over 
all nodes in the networks. For both cases of trapping 
problems, we studied the ATT as an indicator of the 
trapping efficiency. 

We showed that although for a general graph, comput- 
ing ATT to an arbitrary trap is a theoretical challenge, 
for Cayley trees, the ATT to any node can be explicitly 
determined, whose leading term is an increasing function 
of the shortest distance between the central node and the 
trap, implying that the place of trap plays an important 
role in the trapping efficiency. However, for Vicsek frac- 
tals, it is very hard and even impossible to obtain the 
exact expression of ATT for an arbitrary trap, we can 
only obtain the analytic closed-form formula for ATT to 
a peripheral node or the central node, while provided nu- 
merical results for the ATT to other node, with both an- 
alytic and numerical solutions obeying the same leading 
scaling. Thus, different from the case in Cayley trees, the 
trap position has little impact on the trapping efficiency 
for Vicsek fractals. 

For the trapping problem when the trap is uniformly 
distributed, we determined the explicit expressions for 



ATT in both Cayley trees and Vicsek fractals. The ob- 
tained results show that trapping process in Cayley trees 
is much efficient than in Vicsek fractals. In addition, we 
also compared the differences of other aspects for trap- 
ping in the two networks. Finally, we demonstrated that 
the root of all differences for trapping in Cayley trees and 
Vicsek fractals is attributed to their structures. 
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